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It is important to understand the strong external magnetic field generated at the very beginning of 
high energy nuclear collisions. We study the effect of the magnetic field on the charmonium yield and 
anisotropic distribution in Pb+Pb collisions at the LHC energy. The time dependent Schrodinger 
equation is employed to describe the motion of cc pairs. We compare our model prediction of non¬ 
collective anisotropic parameter V 2 of J/tps with CMS data at high transverse momentum. This is 
the first attempt to measure the magnetic field in high energy nuclear collisions. 


It is widely accepted that the most strong magnetic 
field in nature can be generated in the early stage of rel¬ 
ativistic heavy ion collisions. The peaked magnitude of 
the field can reach eB ^ for semi-central Au-|-Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
and 70mJ for semi-central Pb-|-Pb collisions at the Large 
Hadron Collider (LHC) [IM3 , where e is the electron 
charge and m 7 r the pion mass. Considering the interac¬ 
tion between the field and the hot medium formed in the 
later stage of the collisions, one could study fundamen¬ 
tal QCD topological structures 0 in hot and dense nu¬ 
clear matter such as chiral magnetic effect lol- Ha. c hiral 
magnetic wave |13l - [^. chiral separation effect |16l - [lq . 
chiral vortical effect [lOl , chiral electric separation ef¬ 
fect di .and enhancement of elliptic flow of charged par¬ 
ticles mil. Recent experimental results on the study 
of the QCD intrinsic properties at RHIC and LHC can 
be found in Ref. [2|. 

Although we have crude estimation on the strength of 
the initial magnetic field, it may decay very fast and sur¬ 
vive only in the very beginning of heavy ion collisions. 
So far, no experimental determination on the magnetic 
field has been made. In this Letter, we propose using 
the anisotropic production of high momentum charmo- 
nia to probe the initial magnetic field in high energy nu¬ 
clear collisions. Due to their large masses, heavy quarks 
are produced in the early stage through hard scatter¬ 
ings. For low momentum charmonia, they can be pro¬ 
duced in the initial stage and regenerated [ 2 ^-[ 2 | in the 
hot medium, and both will be suppressed through De¬ 
bye screening [3|. High momentum charmonia, on the 
other hand, are purely formed before the formation of the 
hot medium. Except interaction with the magnetic field, 
which will be discussed here, the high momentum char¬ 
monia are almost not affected by the hot medium and 
therefore can be used as an ideal tool to probe the initial 
magnetic field in high energy nuclear collisions. Note, in 
this work we only focus on the spectators induced filed. 
The effects of the delayed decay of the field Q caused by 
the conducting medium will not change our main conclu¬ 
sions. 

Since the life time of the magnetic field ts ^ O.I 
fm/c d is much shorter than the charmorium forma¬ 


tion time t/ ~ 0.5 fm/c from a cc pair to a state 
'k = J,Xc,' ■ the interaction between the cc pair 
and the field could cause some profound effects on (i) 
Charmonium yields: the magnetic field induced force 
will change the charmonium fractions in the cc 

pair |cc) = C^l'k) and thus alters the relative yields 
among different charmonium states, and (ii) Charmo¬ 
nium distribution: the specific direction along the mag¬ 
netic field breaks down the rotational symmetry of the 
system, which leads to an anisotropic charmonium pro¬ 
duction in the transverse plane. Both effects can be ex¬ 
perimentally checked. 

As charm quarks are heavy enough in comparison with 
the inner movement inside a charmonium bound state 
which is determined by the J/ip radius due to the uncer¬ 
tainty principle, rric ~ 1.3 GeV > p ^ l/r = 1/(0.5fm) = 
0.4 GeV, we can ignore the relativistic effect in consid¬ 
ering the inner structure of a charmonium. We employ 
the time dependent Schrodinger equation to describe the 
evolution of a cc pair wave function <I>(t, r^, fg), 

d ^ 

i — <4> = (1) 

dt ^ ' 

with the Hamiltonian operator 


H = 


(Pc - gAg) (pg -t- gAg) 


2mr 


2mr 


V, 


( 2 ) 


where A is the magnetic potential, q = 2e/3 the charm 
quark electron charge, and V the potential between the 
quark and antiquark. We take the Cornell potential to¬ 
gether with the spin-spin interaction [32|, 


V{r) =——+ ar + fie '^’'Sg • Sg. (3) 

r 

By fitting the charmonium spectrum in vacuum, one can 
fix the potential parameters as mg = 1.29 GeV, a = 0.174 
GeV^, a = 0.312, /3 = 1.982 GeV and 7 = 2.06 GeV. 

By separating the cc wave function into a center-of- 
mass part and a relative part $ = with the center- 

of-mass coordinator R = (rg -|- rg)/2 and relative coor¬ 
dinator r = Tg — rg, and further expanding the relative 
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part in terms of the charmonium wave functions 'f'(r), 

1 


$ = 




zPfc-R—2 . ^ 

p 4 mc 




(4) 


the probability |C'^(t)p for the cc pair to be in the char¬ 
monium state 'h satisfies the evolution equation 

f d3r«'*(r)iJs4''(r), (5) 

T , ./ 

\T/f 


where we have separated the Hamiltonian into a vacuum 
part and a magnetic field dependent part, 

H = Ho + Hb, (6) 

with 

2 

AB = -M-B-^(PpxB).r+^(Bxr)2. (7) 

^IiLq ^IIvq 

The charmonium energy and wave function 'l'(r) in 
Eqs. O and ([5]) are determined by Hq, 


Ho'^ = (8) 

Since we have now a special direction, the direction of 
B, the kinetic momentum = P — q{Ac — Ag) with 
the total momentum P = -I- pg is no longer conserved, 

and the conserved momentum is the pseudo momentum 
Pp = P -I- q{Ac — Ag) [ 3 ^. In Eq.([7]) the cc pair interac¬ 
tion with the magnetic field is manifested by three terms: 
the first term is the spin-field interaction with the spin 
magnetic moment fi = q/rricisc ~ Se)j the second term 
comes from the Lorentz force which is proportional to 
the cc momentum, and the third term is the harmonic 
potential which is quadratic in gB and therefore its ef¬ 
fect is much smaller in comparison with the first and 
second terms which are linear in gB. At high momen¬ 
tum, the Lorentz force is the dominant magnetic field ef¬ 
fect. Note that the spin-field interaction makes the spin 
angular momentum no longer a conserved quantity, the 
spin singlet rjc will couple with one of the triplet J/V' 
which carries zero spin component along the magnetic 
field s_B = 0 [H, [13 ■ On the other hand, the Lorentz 
force and the harmonic potential breaks the rotational 
symmetry in the coordinate space. 

Let us take the nuclear colliding direction as the z- 
axis and the impact parameter b parallel to the x-axis. 
While the created magnetic field in heavy ion collisions 
depends on the events, and the magnitude and direction 
of the field fluctuate in space and time EH! , we consider 
here an averaged magnetic field B along the j/—axis in 
the space-time region determined by the colliding energy 
and nuclear geometry, 

{ 2 2 7 ^ 2 :^ 

Bey, 0 < t < t_B, + (b/2)^ < 

0, others. 

( 9 ) 


which, from the relation B = V x A, leads to = Bz 
and Ay = Az = 0. For Pb-|-Pb collisions with centrality 
40% and at LHC energy y/SMN = 2.76 TeV, the geom¬ 
etry parameters and the Lorentz factor are fixed to be 
Ra = 6.6 fm, 6 = 8 fm and 7 c = 1400. From the MC 
simulation Q, we take the magnitude and the life time 
of the magnetic field eB = 25to^ and ts = 0.2 fm/c. At 
RHIC energy, the initially created magnetic field is much 
weaker, eB ~ Q- If the field decays very fast before 
the formation of the hot medium where the Faraday in¬ 
duction may prolong the field Q , the magnetic field effect 
on the charmonium production can be safely neglected. 
In the following we will focus on the magnetic field effect 
at LHC energy. 

To solve the dynamical equation ([5]) for the probability 
magnitude €^,{1), we need the initial value Cii((0) or the 
initial wave function 4>g(0). Suppose the relative motion 
of the cc pair is initially described by a compact Gaussian 
wave package 

$g(0) ~ e "0 (10) 

with averaged relative coordinate rg and distribution 
width (Tq. Since there is no reason to choose a special 
direction before the magnetic field is introduced, the az¬ 
imuth angles of rg = rg(sin0g cos ()ig, sin0g sin(/)g, cos 0g) 
are randomly distributed, and we do ensemble average 
over 9q and (j)Q in the final state calculations. The re¬ 
maining two parameters rg and erg can be determined by 
fitting the charmonium fractions in p-|-p collisions. Sup¬ 
pose an initial point-like wave function 5(r — r') develops 
as ) /C ‘ ) with a constant expansion velocity v, the 
initial Gaussian wave package evolves as 

f (r-r')2 

$r.(t) ^ "0 e ^ e (11) 


in p-|-p collisions. It is always a Gaussian wave pack¬ 
age but with a time dependent width a^=aQ+ v^t^. 
Calculating the probability |C' 5 -(t/)p = |('I'|$r(t/))p in 
vacuum with H = 0 at the charmonium formation time 
tf, and taking the ex per imentally observed decay branch 
ratios H(4' —>• J/ip) we obtain the fractions of the 
direct J/ip production and feed down from the excited 
states in the finally observed prompt J/ips, 


R\jj{t) 


|C^(t)|2H(^ ^ J/^) 

^ j/i;) 


( 12 ) 


with the definition of B{J/'ip —>■ J/ip) = 1- Taking 
R{J/iIj) = 60%, i?(r/>') = 10% and R{^ = 30% from 
the recent p-l-p data at LHC energy |36l - l3^ and the char¬ 
monium formation time tf = 0.5 fm/c [3ll |. the remain¬ 
ing two parameters in the initial wave package $,.(0) are 
fixed to be rg = 0.68 fm and erg = 0.02 fm, which corre¬ 
spond to the expansion velocity v = 0.72c and the width 
atf = 0.38 fm. 

With the known initial wave function (0) or the ini¬ 
tial probability magnitudes Cj/^{0) : C'^'(O) : C^^{0) = 
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1 : 0.5 : 1.2, we can solve the evolution equation ([5]) for 
Cq,{t) in the time region 0 < t < ts when the magnetic 
field exists. After that the wave function evolves accord¬ 
ing to the expansion rule (ITT]) with r) as the initial 

condition, 


for t>tB- 

(13) 

We now show our numerical results about the mag¬ 
netic field effect on the charmonium production in heavy 
ion collisions. We focus on the central rapidity region in 
Pb-|-Pb collisions with impact parameter 6 = 8 fm and 
at LHC energy ^Jsnn = 2.76 TeV. In our calculation, 
the initial cc pairs which are determined by the colliding 
energy and the nuclear geometry are assumed to be dis¬ 
tributed homogeneously in the overlapped region of the 
two colliding nuclei. To have a good precision in solving 
the coupled equations ([5]) for C'ii.(t), we take a cutoff in 
the sum over the charmonium eigenstate 'h: the main 
quantum number n < 6 and the orbital angular momen¬ 
tum number I < 7. 


Pb+Pb Collisions at \/Sfj^| = 2.76 TeV 



FIG. 1: (Color online) The time evolution of J/ijjs from di¬ 
rect production (solid lines) and feed down from ip' (dashed 
lines) and Xc (dot-dashed lines). The thick and thin lines 
represent the results from the calculations with and without 
the external magnetic filed, respectively. As indicated by the 
vertical short-dashed line, the magnetic field only lasts during 
the time t <tB ~ 0.2 fm/c. 


In FiglU the J/ip fractions from different chan¬ 

nels are shown as functions of cc evolution time. The 
results with and without the external magnetic field are 
displayed by thick and thin lines, respectively. These 
fractions are extracted at fixed transverse momentum 
Pt = 10 GeV/c from Pb-|-Pb collisions at impact pa¬ 
rameter 5 = 8 fm and LHC energy v^^vAf = 2.76 TeV. 
In case there is no magnetic field, the direct J/ip produc¬ 
tion and feed down from ip' are slightly enhanced, while 
the contribution from Xc is weakly suppressed, indicat¬ 
ing a Xc decay to J/ip and ip' with quantum selection 
rule 61 = \. When the magnetic field is turned on, the 
converting from the high lying state Xc to J/ip and ip' 
becomes much more dramatic, as one can see the thick 
lines in FigHJ Since the magnetic field introduces a spe¬ 
cial direction, the cc motion becomes anisotropic. Here 
we fix the cc azimuthal angle in the transverse plane 
ip = arctan(py/p 2 :) = 0 where the magnetic field ef¬ 
fect is expected to be the strongest. In the time period 
t < tBi the external magnetic field serves as a stimulator 
that enhances the quantum mechanics allowed transition 
from Xc to J/ip and ip'. After the field is off at t > ts, 
the variation of all the fractions with time becomes mild 
again. At the formation time t/ = 0.5 fm/c, the relative 
enhancement for both direct J/ip production and feed 
down from ip' are found to be 10%, and the contribution 
from Xc decay is relatively suppressed by 23%. 

Having discussed the different contributions to J/ip 
production, we now look at the yields for J/ip, ip' and 
Xc- For J/ip, the ratio between the yields and 
with and without the magnetic field can be expressed in 
terms of the probabilities and \C^\tf)\'^ and 


Pb-t-Pb Collisions at \/sT| = 2.76 TeV 
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FIG. 2: (Color online) The magnetic field induced anisotropic 
behavior of the fractions A® in the J/ip yield at the formation 
time tj. 
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(14) 


Using the above calculated probabilities and fractions 
at LHC energy, the maximum J/ip enhancement at 
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= 0 and pt = 10 GeV/c is 13%. For the excited 
states, neglecting the feed down from the higher eigen 
states of 7^0, we have the yield ratios /N^) = 
= 1.29 for and = 

l^'xf= 0.84 for Xc, which mean, due 
to the magnetic held, a relative enhancement of 29% and 
suppression of 16% for ■0' and Xc productions in Pb+Pb 
collisions at the LHC energy. 

Fig 13 shows the magnetic held induced anisotropic 
property of the fractions Rq,{tf) in the transverse plane 
at charmonium formation time tj. The strength of the 
Lorentz force acting on the charmonia is most strong at 
(p = 0, and decreases monotonously with increasing az¬ 
imuthal angle ip. Finally at p = 1^12, the force disap¬ 
pears and only the weak harmonic potential exists, the 
fractions approach to their vacuum values. 


Pb+Pb Collisions at = 2.76 TeV 



Transverse Momentum (GeV/c) 


Pb+Pb Collisions at \/S|0 = 2.76 TeV 



FIG. 3: (Color online) The magnetic field induced transverse 
momentnm dependence of the fractions in the J/0 yield 
at the formation time t/. 


FIG. 4: (Color online) The transverse momentum depen¬ 
dence of J/0 V 2 - The solid and dashed lines are the initially 
produced non-collective J/0 V 2 with and without magnetic 
field, and the data at high pt are from the CMS collabo¬ 
ration [^. As a comparison, the collective J/0 V 2 from a 
transport model is shown as the dot-dashed line. 


lective motion of the medium created in high energy nu¬ 
clear collisions [39l - l4ll| . At the LHC energy, J/0 V 2 has 
been reported in the region of pr < 10 GeV/c • Since 
the collective motion is dominated by the bulk interac¬ 
tions in relatively low momentum region, those high pp 
charmonia generated in the initial stage are not expected 
to be sensitive to the nature of the hot medium. How¬ 
ever, the Lorentz force induced anisotropic production in 
the transverse plane, shown in FiglJl may result in a non¬ 
collective J/0 V 2 at high Pt- The J/0 V 2 as a function 
of rapidity and transverse momentum is defined as 


We now show in Fig|3]the magnetic field induced trans¬ 
verse momentum dependence of the fractions at 

charmonium formation time tf. We again fix the angle 
(/? = 0 to see the maximum magnetic field effect. At 
Pt = 0 the Lorentz force disappears, the slight deviation 
of the fractions from the vacuum values comes from the 
weak harmonic potential. The strength of the Lorentz 
force increases with pp, and the change of the fractions 
for all channels becomes larger till pp ^ 3 GeV/c. In 
higher pp region, the change stops and the transverse 
momentum dependence seems flat. Although the force 
is increasing with pp, the higher the pp the shorter the 
time the cc pair experiences the field. The combination 
of the momentum dependence of the Lorentz force and 
the finite time of the external magnetic field leads to the 
observed saturation in Figl31 

Usually, the observed event anisotropy V 2 is discussed 
in the framework of hydrodynamics, representing the col- 


V2{V,Pt) 


Jo^ Nj/^{r],pt, p) cos{2p)dp 
Jo^ Nj/,^{r],pt,p)dp 


(15) 


The numerical result at central rapidity is shown in Fig 0] 
In normal calculations without considering the magnetic 
field [39l - l4ll| . the initial J/0s are isotropically produced 
with t ;2 = 0 (dashed line). However, due to the mag¬ 
netic field induced Lorentz force, the high pp J/ijjs ac¬ 
quire sizable non-collective V 2 (solid line), which explains 
reasonably well the GMS data for prompt J/ijjs in 
the high pp region. As a comparison. We show in Figl4] 
also the collective V 2 (dot-dashed line) calculated from a 
dynamical transport model [s^. Different from the col¬ 
lective flow which comes from the J/'ijj regeneration at 
low and intermediate pp, the non-collective V 2 induced 
by the magnetic field is mainly in high pp region. As one 
can see in Figs[Tl[3and[31 Xc is suppressed, corresponding 
to the J/0 and 0' enhancement. Hence it is interesting 
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to point out that if the observed non-collective V 2 for J/'0 
is positive, the high px for Xc should be negative. 

In summary, by solving the time dependent 
Schrodinger equation for the cc pairs, we discussed 
the effects of the initially created magnetic field on 
the relative yields and anisotropic distributions of the 
charmonium states in high energy nuclear collisions. 
For Pb-|-Pb collisions at LHC, we found (i) the directly 
produced and ip’ decayed J/ips, are enhanced, while 
the J/ips from Xc decay are suppressed, and (ii) the 
non-collective J/ip V 2 at large transverse momentum. 


Pt > 8 GeV/c, is as large as 4% and comparable with 
the CMS data. We wish to point out that in high energy 
nuclear collisions the initial magnetic field is related to 
several important measurements reflecting the intrinsic 
structure of the QCD. The measurement of the strong 
J/ip anisotropy at large px can be used to diagnose the 
magnetic filed in such collisions. 
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